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Abstract The angular momentum of a star is an important astrophysical quan- 
£Nlj . tity related to its internal structure, formation, and evolution. Helioseismology 

yields S Q = 1.92 x 10 41 kg m 2 s" 1 for the angular momentum of the Sun. We 
show how it should be possible to constrain it in a near future by using 
the gravitomagnetic Lense-Thirring effect predicted by General Relativity for 
the orbit of a test particle moving around a central rotating body. We also 
discuss the present-day situation in view of the latest determinations of the 
supplementary perihelion precession Acb^ of Mercury A fit by Fienga et al. 
(Celestial Mech. Dynamical Astron. Ill, 363, 2011) of the dynamical models of 
several standard forces acting on the planets of the solar system to a long data 
q-( record yielded Ari/™ eas ' = 0.4 ± 0.6 milliarcseconds per century. The modeled 

$h . forces did not include the Lense-Thirring effect itself, which is expected to be 

as large as cbj^ = -2.0 mas cty -1 from helioseismology-based values of S s . 

By assuming the validity of General Relativity, from its theoretical prediction 
for the gravitomagnetic perihelion precession of Mercury one can straightfor- 
wardly infer S Q < 0.95 X 10 41 kg m 2 s -1 . It disagrees with the currently available 
values from helioseismology. Possible sources for the present discrepancy are 
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^■j. ■ examined. Given the current level of accuracy in the Mercury ephemerides, 

the gravitomagnetic force of the Sun should be included in their force models. 
MESSENGER, in orbit around Mercury since March 2011, will collect science 
data until 2013, while BepiColombo, to be launched in 2015, should reach Mercury 
in 2022 for a year-long science phase: the analysis of their data will be important 
in effectively constraining S Q in about a decade or, perhaps, even less. 
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1. Introduction 

The angular momentum [S] of a main-sequence star can yield relevant informa- 
tion about its inner properties and its activity (Antia, Chitre, and Thompson, 
2000; Yang and Bi, 2006, 2008; Bi et al, 2011). More specifically, S is related to 
the stellar internal rotation rate (Pijpers, 1998). Moreover, it is an important 
diagnostic for testing theories of stellar formation (Pijpers, 2003). The angular 
momentum can also play a decisive role in stellar evolution, in particular to- 
wards higher mass (Tarafdar and Vardya, 1971; Wolff, Edwards, and Preston, 
1982; Vigneron et al, 1990; Wolff and Simon, 1997; Di Mauro, Murabito, and 
Paterno, 2000; Komm et al, 2003; Herbst and Mundt, 2005; Jackson, MacGregor, 
and Skumanich, 2005; Antia, Chitre, and Gough, 2008). 

In this work, we propose to put dynamical constraints on S , independent of 
any model of the Sun's interior, by exploiting General Relativity. See also Iorio 
(2008) for earlier investigations. The plan of the article is as follows: In Section 2 
we summarize the current level of knowledge of the solar angular momentum 
from other, non-dynamical techniques. In Section 3 we introduce the general- 
relativistic Lense-Thirring effect, which is related to the angular momentum. 
We apply it in Section 4 in connection with the latest determinations of the 
orbital motions of the inner planets of the solar system (Fienga et al, 2011). In 
Section 5 we discuss the impact that other competing dynamical effects may 
have. Section 6 summarizes our findings. 

2. The Angular Momentum of the Sun: Present-Day Knowledge 

The angular momentum is one of some global properties of the Sun of as- 
trophysical interest, which are related to the internal rotation rate through 
integral equations (Pijpers, 1998). The solar internal rotation rate can be de- 
termined using inverse techniques of helioseismology (Christensen-Dalsgaard, 
2002; Kosovichev, 2003; Di Mauro, 2003, 2008) applied to data collected from 
both the ground-based Global Oscillation Network Group (GONG: Harvey et ah, 
1996) and the instruments carried onboard the SOlar and Heliospheric Observa- 
tory (SOHO: Domingo, Fleck, and Poland, 1995), such as the Solar Oscillations 
Investigation/Michelson Doppler Imager (SOI/MDI: Scherrer et al, 1995) and Global 
Oscillations at Low Frequencies (GOLF: Gabriel et al, 1995), and the Solar Dynamics 
Observatory/Helioseismic and Magnetic Imager (SDO/HMI: Scherrer et al, 2012). 

Such observations pertain the solar oscillation frequencies. The observed 
solar oscillations correspond to standing acoustic waves maintained by pressure 
forces, which form the class of the p-modes, and to standing surface gravity 
waves, maintained by gravity, known as /-modes. In addition, one should 
mention the probable, although quite discussed, existence of the internal gravity 
waves, g-modes, which are sensitive to the structure and rotation of the deeper 
interior of the Sun (Di Mauro, 2008). Oscillations have several advantages over 
all the other observables: their frequencies can be measured with high accuracy, 
and depend in a quite simple way on the equilibrium structure of the Sun (Di 
Mauro, 2008). 
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Table 1. Values of the solar angular momentum, in units of 10 41 kg m 2 s _1 , retrieved from 
the literature. The uncertainties, when released by the authors, have been quoted as well. 
With the exception of the figure by Livingston (2000), inferred from surface rotation, all 
the other values quoted come from helioseismology. Pijpers (1998) used data from both 
GONG and SOHO/MDI. Di Mauro, Murabito, and Paterno (2000) relied upon MDI and 
GOLF. Komm et al (2003) used GONG and MDI data as well, but in a different way with 
respect to Pijpers (1998) and Antia, Chitre, and Thompson (2000). The values by Yang and 
Bi (2006, 2008) and Bi et al (2011) were obtained from models taking into account the 
effects of the magnetic field: we quoted only those figures consistent at 1 - 3a level with 
helioseismological results. 



S [10 41 kgn^s" 1 ] 


a Se [rO^kgi^s" 1 ] 


Reference 


1.90 


0.015 


(Pijpers, 1998) 


1.63 




(Livingston, 2000) 


2.02 


0.04 


(Di Mauro, Murabito, and Paterno, 2000) 


1.91 




(Antia, Chitre, and Thompson, 2000) 


1.94 


0.05 


(Komm et al, 2003) 


1.91 




(Yang and Bi, 2006) 


2.045 




(Yang and Bi, 2008) 


2.02 




(Bi et a?., 2011) 



Table 1 summarizes the present-day knowledge of S Q according to mainly 
helioseismology. On average, it is 

S = 1.92xl0 41 kgm 2 s _1 . (1) 



3. Using General Relativity: the Lense-Thirring Effect 

Contrary to the Newtonian theory of gravitation, the Einsteinian General Rel- 
ativity predicts the existence of dynamical effects which can, in principle, be 
used to measure, or, at least, constrain the angular momentum of a rotating 
body since their analytical expressions are directly proportional to S. 

The exterior spacetime metric of a typical astronomical object can adequately 
be described within the Parameterized Post-Newtonian (PPN) approximation 
(Soffel, 1989). Within such a framework, a stationary gravitomagnetic field arises 
around a slowly rotating mass with proper angular momentum S (Thorne, 
MacDonald, and Price, 1986; Thorne, 1988; Mashhoon, Gronwald, and Licht- 
enegger, 2001). Gravitomagnetism does not refer to any combined effect of 
electromagnetism and gravitation. Such a purely formal denomination is due to 
the fact that the gravitomagnetic field affects a test particle moving with velocity 
v with a noncentral, Lorentz-like acceleration (McCarthy and Petit, 2004) which 
is analogous to the one felt by a moving electric charge in a magnetic field in 
the framework of Maxwellian electromagnetism. It is the cause of the so-called 
Lense-Thirring effect (Lense and Thirring, 1918) consisting of small secular 
precession of some of the osculating Keplerian orbital elements of the orbit of a 
test particle. 
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The gravitomagnetic precession of a gyroscope (Pugh, 1959; Schiff, 1960) was 
recently measured in a dedicated space-based experiment, known as Gravity 
Probe B (GP-B: Everitt, 1974), performed with four cryogenic gyroscopes carried 
onboard a drag-free spacecraft orbiting the Earth. The prediction of General 
Relativity for such an effect was successfully corroborated at a claimed 19% 
accuracy (Everitt et ah, 2011). It may be interesting to recall that Haas and Ross 
(1975) proposed to measure the same effect with a drag-free spacecraft orbiting 
the Sun. 

In this article, we will show how to use the gravitomagnetic orbital pre- 
cessions to constrain the solar angular momentum. Some spacecraft-based 
missions were proposed in the more or recent past to accurately measure the 
Sun's gravitomagnetic field by means of its direct effects on the propagation of 
electromagnetic waves. They are the Laser Astrometric Test of Relativity (LATOR: 
Turyshev, Shao, and Nordtvedt, 2004), which aimed to directly measure the 
frame-dragging effect on the light with a » 0.1% accuracy (Turyshev et ah, 2009), 
and the Astrodynamical Space Test of Relativity using Optical Devices I (ASTROD I: 
Ni, 2008), whose goal is to measure the gravitomagnetic component of the time 
delay with a 10% accuracy (Appourchaux et al, 2009). We notice that any tests 
of general-relativistic solar frame-dragging would necessarily be accurate to a 
percent level. Indeed, the angular momentum of the Sun should be considered 
as known independently of the Lense-Thirring effect itself; helioseismology 
yields results accurate just to that level, as shown by Table 1. 

Iorio (2012a) analytically worked out the secular Lense-Thirring precessions 
of the Keplerian orbital elements for a generic orientation of the angular mo- 
mentum of the central body in a form that allows for easy comparisons with the 
observation-based quantities usually determined by the astronomers. They are 
(Iorio, 2012a) 

0, 
0, 

2GS(k x cos Q+fc y sin O) 
c 2 a 3 (l-e 2 ) 3 ' 2 ' 

(2) 

2GS[L+cot/(£ y cosn-^ sin £2)] 
c 2 a 3 (l-e 2 ) 3 ! 2 ' 

2GS{2[£ 2 cos Z+sin l(k x sinQ-^cos Q)]-[L sinZ+cosl(fc v cosQ-£ r sinO)] tan(Z/2)| 

c 2 a 3 (l-e 2 ) 3 ' 2 ' 

0; 

(k x , k y , is the unit vector of the rotation axis of the central body 
in the reference frame adopted, while a, e, I, Q, tD, and M are the semimajor 
axis, the eccentricity, the inclination, the longitude of the ascending node, the 
longitude of pericenter, and the mean anomaly, respectively, of the orbit of the 
test particle (Murray and Dermott, 1999). In particular, I is the inclination of the 
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orbital plane to the reference (x, y) plane, Q is an angle in the reference (x, y) 
plane counted from the reference x-direction to the line of the nodes which is 
the intersection of the orbital plane to the reference (x, y) plane, and cD = Q + o> 
is a "dogleg" angle since the argument of pericenter co is an angle in the orbital 
plane counted from the line of the nodes to the point of closest approach (Murray 
and Dermott, 1999). The precessions of Equation (2) are averages over one full 
orbital revolution of the test particle; they were obtained (Iorio, 2012a) by using 
the Lagrange perturbative equations (Murray and Dermott, 1999) applied to a 
suitable perturbing potential (Barker and O'Connell, 1975) averaged over one 
orbital revolution. It should be noticed that almost all of the usual derivations 
of the Lense-Thirring effect existing in literature are based on the particular 
choice of aligning S along the z-axis. For some exceptions, in different contexts, 
see Barker and O'Connell (1970) and Will (2008). However, the authors of such 
works either did not explicitly work out the precessions or adopted a different 
parameterization in such a way that their results could not straightforwardly 
be adapted to the present context. 

In the case of the Sun, Fienga et al. (2011) adopted the mean Earth's Equator 
at the epoch J2000.0 for their analysis; in such a frame, the unit vector of the 
Sun's spin axis [k Q ], entering Equation (2), is not aligned with the z-axis. Indeed, 
it turns out 

fc° = cos 6 cos a Q = 0.122, 

= cos 6© sin a Q = -0.423, (3) 

kf = sin do = 0.897, 
where (Seidelmann et al., 2007) 

<5 = 63.87 deg, 

(4) 

a Q = 286.13 deg, 

are the declination and the right ascension, respectively, of the Sun's North Pole 
of rotation with respect to the mean terrestrial equator at J2000.0. 

4. Constraints from Planetary Orbital Motions 

Recently, Fienga et al. (2011) processed a huge amount of planetary observations 
of various kinds, covering almost one century (1914-2010), with the dynamical 
force models of the INPOPlOa ephemerides (Fienga et at, 2010). The data in- 
cluded, among other things, also high-quality Doppler range-rate observations 
to the MErcury Surface, Space ENvironment, GEochemistry, and Ranging spacecraft 
(MESSENGER: Solomon et al, 2007) collected 1 during its three flybys of Mercury 



1 More precisely, Fienga et al. (2011) neither used direct observations to MESSENGER to fit the IN- 
POPlOa ephemerides nor estimated the spacraft's orbit. Instead, they used normal points extracted 
from the SPICE/NAIF MESSENGER orbit [http://naif.jpl.nasa.gov/naif/]. Fienga et al. (2011) plan to 
themselves analyze such data in the coming months. 
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in 2008-2009 (Smith et at, 2010). Thus, our knowledge of the orbit of the 

innermost planet of the solar system has been greatly improved. Generally 

speaking, Mercury is a very difficult planet to observe from the Earth because 

of a number of reasons (Balogh et at, 2007). Also orbiting Mercury by spacecraft 

is a challenge, since the planet is deep inside the gravitational potential well 

of the Sun (Balogh et at, 2007). Moreover, a very hostile thermal environment 

awaits any spacecraft in Hermean orbits (Balogh et at, 2007). As a result, to 

date only one spacecraft, Mariner 10, reached Mercury more than 30 years ago. 

See Balogh et al. (2007) for a general overview on the missions to Mercury. 

MESSENGER was inserted into orbit around Mercury on 1 7 March 2011 (Harris, 

2011) for a one-year-long, near-polar-orbital observational campaign (Solomon 

et ah, 2007). A year-long extension until March 2013 was recently approved 

[http://www.space.com/13655-mercury-spacecraft-messenger-mission-extension.html].! 

BepiColombo (Grard and Balogh, 2001; Benkhoff et at, 2010) is another planned 

mission to Mercury whose launch is currently planned for 2015 [http://sci.esa. int/science-e/www/area/index.c 
while its one-year mission around Mercury should start after its arrival sched- 
uled for 2022. Within the framework of the Mercury Orbiter Radioscience Experi- 
ment (MORE: less, Asmar, and Tortora, 2009), accurate measurement of its range 
should be able to notably improve the determination of Mercury's orbit as well 
(less and Asmar, 2007). 

Apart from usual Newtonian mechanics, including the Sun's quadrupole 
mass moment J2 (Rozelot and Damiani, 2011) as well, the mathematical mod- 
els of INPOPlOa for the solar system dynamics included also the first post- 
Newtonian (1PN) static, Schwarzschild-like component of the gravitational 
field of the Sun (McCarthy and Petit, 2004) expressed in terms of the usual 
Parameterized-Post-Newtonian (PPN) parameters y and ji, which are equal 
to one in General Relativity. Instead, Fienga et at (2011) did not model the 
gravitomagnetic field of the Sun. 

Fienga et al. (2011) estimated, among other things, corrections AQ (meas ' and 
^^(meas) j- fag standard Newtonian-Einsteinian secular precessions of the lon- 
gitudes of node Q and perihelion ffl for the first six planets. In principle, 
such corrections account for any unmodeled/mismodeled dynamical features 
of motion, so that they could be used to pose constraints on the magnitudes of 
putative anomalous effects, if any, with respect to the modeled ones. Concerning 
Mercury, they are (Fienga et at, 2011) 

A ^(meas) _ 1-4 + 1-g mas cty -l 

(5) 

Ad) (meas) _ .4 ± 0.6 mas cry" 1 , 

where mas cty -1 is a shorthand for milliarcseconds per century. The values of 
Equation (5) are statistically compatible with zero, so that they could be used, 
in principle, to put upper bounds on S by comparing Equation (2) to Equation 
(5). It makes sense since the predicted Lense-Thirring precession for it is larger 
than 0.6 mas cty" 1 , as it turns out by using the values for S listed in Table 
1. In Table 2 we display the corresponding Lense-Thirring precessions for the 
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Table 2. Expected Lense-Thirring perihelion precession for Mercury, in mas cty 1 , ac- 
cording to Equation (2) computed with the figures of Table 1 for the solar angular momentum. 
See also (Iorio, 2008). They are larger than the uncertainty released by Fienga et a]. (2011) for the 
correction to the standard perihelion precession of Mercury, which amounts to 0.6 mas cty -1 . We 
also show the discrepancy among the predicted and the determined Affljj 1638 ' of Equation 
(5) in number of a. 



ffl^ T ' (mas cty 1 ) 


Discrepancy (number of a) 


Reference for S Q 


-2.0 


4 


(Pijpers, 1998) 


-1.7 


3.5 


(Livingston, 2000) 


-2.1 


4.2 


(Di Mauro, Murabito, and Patemo, 2000) 


-2.0 


4 


(Antia, Chitre, and Thompson, 2000) 


-2.0 


4 


(Komm et ah, 2003) 


-2.0 


4 


(Yang and Bi, 2006) 


-2.1 


4.2 


(Yang and Bi, 2008) 


-2.1 


4.2 


(Bi et ah, 2011) 



perihelion of Mercury. On average, they yield (Iorio, 2008) 



ci£ T) = -2.0 mas cty" 1 ; (6) 



departures from such a value are not relevant since they are smaller than 0.6 
mas cty -1 . Iorio (2008) used the usual formula with S parallel to the z-axis, but 
it turns out that, contrary to the node and the inclination, it does not affect the 
perihelion precession to a significant level of accuracy. It can be noticed that there 
is a discrepancy of 4a, on average, among the predictions of General Relativity, 
computed using the helioseismically-inferred values for S Q , and the admissible 
range -0.2 mas cty" 1 < A(i> ( ™ eas ' < lmas cty" 1 for unmodelled dynamical effects 

affecting the perihelion of Mercury obtained from observations by Fienga et al. 
(2011). 

We did not consider the Lense-Thirring node precessions since they are 

expected to be smaller than 1 mas cty" 1 . 

A direct comparison of Equation (2) to Equation (5) yields, in units of 
10 41 kg m 2 s" 1 , 

S < 0.95. (7) 

It is notably smaller that the figures of Table 1, at more than 10cr level. 
In Section 5 we offer a preliminary discussion of possible explanations. 
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5. Possible Explanations 

It must be stressed that the previous considerations would be valid if the sup- 
plementary advance of the perihelion of Mercury of Equation (5) were entirely 
explained in terms of frame-dragging. 

Actually, in principle, there is the possibility that the unmodeled gravito- 
magnetic effect was partially or totally removed from the post-fit signature in 
the data-reduction process, having been somewhat "absorbed" in the values of 
some of the standard parameters estimated in the fits like, e.g., the planetary 
initial conditions. After all, the expected magnitude of the Lense-Thirring effect 
is about of the same order of magnitude of the present-day accuracy in deter- 
mining the orbit of Mercury. The PPN parameters f> and y entering the modeled 
static 1PN solar field were kept fixed to their general-relativistic values in the 
solutions yielding A(i>< meas > and AQ< meas > (Fienga et al, 2011). 

A possible conventional explanation in terms of orbital dynamics might be 
the following. The corrections Acb were estimated by Fienga et al. (2011) by 
modeling, among other things, the Newtonian action of the quadrupole mass 
moment [J2] of the Sun (Rozelot and Damiani, 2011) as well: more specifically, 

was kept fixed to a reference value, which seems to be (Fienga et al., 2010) 

J® = 1.8 X 10" 7 . (8) 

For a historical overview of the impact of J2 on the measurement of the larger 
1PN Schwarzschild-like perihelion precession of Mercury (e.g., Pireaux and 
Rozelot, 2003). Since such a physical quantity is known with a necessarily limited 
accuracy, of the order of (Fienga et at, 2010, 2011) 10%, the corrections Acb 
account, in principle, not only for the completely unmodeled Lense-Thirring 
effect, but also for the mismodeled precessions Acb'? due to the solar oblateness 
itself. Thus, a mutual cancelation might have occurred leaving just Equation (5). 
To check this possibility, one has to consider the sum of the precessions of the 
perihelion of Mercury caused by the Lense-Thirring effect and by a correction 
A /? to the reference value used; analytical expressions for the long-term orbital 
precessions induced by the oblateness of the primary for an arbitrary direction 
of k were recently obtained in a form suitable for comparison with observations 
in terms of usual osculating Keplerian orbital elements (Iorio, 2011). 

Similar considerations hold also for the 1PN Schwarzschild-like component 
of the Sun's gravitational field: in the PPN formalism, the resulting perihelion 
precession for Mercury consists of the familiar 43.98 arcsec cty -1 rescaled by 
the multiplicative coefficient (2 + 2y - jS)/3 written in terms of the usual PPN 
parameters |8 and y. Actually Fienga et al. (2011) modeled it in the solution 
yielding Equation (5) by keeping |8 and y fixed to one, but the current uncertainty 
in f>, evaluated by Fienga et al. (2011) themselves as large as 

jS - 1 = (-0.41 + 0.78) x 10" 4 (9) 

in a different solution by using the bounds on y from Bertotti, less, and Tortora 
(2003), may translate into a further mismodeled competing perihelion preces- 

sion Act)' which, added to Ad)* and the expected Gyrf , may yield just Equation 
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(5) from a mutual cancelation. It turns out that A©^ alone, computed with 
Equation (9) and added to cb^ T , is not able to reproduce Equation (5). Instead, 

JO 

by including also Atb^, computed from Equation (20) - Equation (21) of Iorio 
(2011) with 

0.1 x 10 -7 < A/f < 0.2 x 10 -7 , (10) 

a range compatible with Equation (5) is obtained. 

In principle, also the impact of the oblateness of Mercury itself on its orbital 
motion should be considered, since it was not modeled in INPOPlOa. However, 
it turns out that it is insufficient to cancel the Lense-Thirring effect to the level 
of Equation (5). Indeed, by using Equation (20)-Equation (21) of Iorio (2011) 
with (Seidelmann et al, 2007) 

= 2.4397 xl0 6 m, 

a H = 281.01 dee, 
? h (11) 

5y = 61.45 deg, 



and (Smith et al, 2010) 



jf = (1.92 ±0.67) xlO -5 , (12) 



the resulting perihelion precession is just 

cb 7 * = (0.03 ±0.01) mas cry" 1 . (13) 

In principle, viable candidates are also the major asteroids and the ring of 
the minor asteroids, although their dynamical action was accurately modeled 
by Fienga et al. (2011). Indeed, both individual distant bodies and an external 
massive ring induce secular perihelion precessions. Let us, first, consider the 
impact of the uncertainty in the mass of some of the individual major asteroids 
such as (1) Ceres, (2) Pallas and (4) Vesta. They are at about 2.4 — 2.8 au. The 
present-day uncertainty in their masses is (Luzum et al., 2011) 

cr fflast = 3 x 1O" 12 M ; (14) 

Fienga et al. (2011) report for them similar or smaller uncertainties. Thus, 
Equation (13) -Equation (14) of Iorio (2012b) yield 

AtbS* * 1 x 10 -3 mas cry" 1 , (15) 

which is, actually, of no concern. The secular perihelion precession induced by 
an external massive ring with mass m T and radius R r has been treated by a 
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number of authors with different approaches (e.g., Fienga et ah, 2008; Kuchynka 
et ah, 2010; Iorio, 2012c). For (Kuchynka et al, 2009) 

R r = 2.80 AU, o„ h = 2 x 10- n M O/ (16) 

they yield 

Atbj^ w 0.02 mas cty -1 , (17) 

which is negligible as well since it is 30 times smaller that the accuracy in 
determining the supplementary perihelion advance of Mercury. 

Finally it is not unreasonable to argue that, ultimately, Equation (5) is based 
just on three normal points from MESSENGER, so that the entire matter should 
be left on hold until more data from the current Mercury orbiter will be gathered 
and analyzed. 

6. Summary and Conclusions 

We showed how it is possible, in principle, put constraints on the angular 
momentum of the Sun, independently of models of its interior, by exploiting the 
general-relativistic gravitomagnetic Lense-Thirring effect. It consists of secular 
precessions of the node and the pericenter of a test particle orbiting a central 
rotating body which are proportional just to its angular momentum. 

Contrary to the 1PN gravitoelectric Schwarzschild-like component of the 
gravitational field of the Sun and all standard Newtonian effects, the solar 
gravitomagnetic field was not included in the mathematical models of the forces 
acting on the planets of the solar system recently fit to long data series by 
a team of astronomers. From the observations, including also three flybys of 
MESSENGER of Mercury, they computed the maximum allowed range of values 
for any unmodeled/mismodeled effect impacting the secular precessions of the 
planetary perihelia. 

If one used them to straightforwardly infer upper bounds on the Sun's angu- 
lar momentum from a comparison with the theoretical prediction for the Lense- 
Thirring perihelion precession, it would turn out S < 0.95 x 10 41 kg m 2 s . 
Such a result would be in nice disagreement with the non-dynamical val- 
ues for S e inferred from helioseismology, which, on average, point toward 
S® = 1.92 X 10 41 kg m 2 s -1 . We discussed some possible explanations. In particu- 
lar, we looked at several competing dynamical effects that may have conspired 
to such an outcome. A viable candidate is the residual precession due to the 
mismodeling in the solar quadrupole mass moment ]i, currently known with 
a « 10% accuracy. Other potential sources of aliasing such as the oblateness of 
Mercury itself, the major asteroids, and the ring of minor asteroids are ruled out 
because their effects are smaller than the current accuracy in determining the 
orbit of Mercury. 

Since the Lense-Thirring perihelion precession of Mercury, as expected from 
General Relativity with helioseismology data, is -2.0 mas cty -1 and the present- 
day accuracy in constraining the secular perihelion precession of Mercury is 
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0.6 mas cty -1 from INPOPlOa ephemerides, it seems possible to effectively 
constrain the solar angular momentum in the near future. To this aim, the 
Sun's gravitomagnetic force should be explicitly included in the software used 
by astronomers to reduce data by modifying the coded dynamical force models; 
it should not be a prohibitive, time-consuming task. As a first step, an analysis of 
the complete data record from MESSENGER, which was inserted in orbit around 
Mercury in March 2011 for a year-long science phase extended until March 
2013, could be implemented in the next few years. Over a longer timescale, 
BepiColombo, whose launch from the Earth and arrival to Mercury are scheduled 
for 2015 and 2022, respectively, will further improve our knowledge of the orbit 
of Mercury allowing for a more refined analysis. It is hoped that astronomers 
engaged in the production of modern ephemerides will undertake such an 
effort, which may be rewarding in no more than ten years. 
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